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Abstract 

Two-dimensional random tilings of rhombi can be 
seen as projections of two-dimensional membranes 
embedded in hypercubic lattices of higher dimen- 
sional spaces. Here, we consider tilings projected 
from a D-dimensional space. We study the limiting 
case, when the quantity D, and therefore the num- 
ber of different species of tiles, become large. We 
had previously demonstrated that, in this limit, 
the thermodynamic properties of the tiling become 
independent of the boundary conditions. The exact 
value of the limiting entropy and finite D corrections 
remain open questions. Here, we develop a mean- 
field theory, which uses an iterative description of 
the tilings based on an analogy with avoiding ori- 
ented walks on a random tiling. We compare the 
quantities so-obtained with numerical calculations. 
We also discuss the role of spatial correlations. 



1. Introduction 

The understanding of the stability of quasicrystalline 
materials, discovered about 15 years ago moti- 
vates the study of random tilings. Indeed, it is be- 
lieved that in these materials, atomic motifs form 
geometric tiles; the random tiling hypothesis states 
that these tiles can be rearranged via local degrees 
of freedom, which gives access to a finite config- 
urational entropy which is supposed to favor the 



quasicrystalline symmetry against other competitive 
crystalline phases Q. Indeed, such symmetries, a 
priori inherent to Penrose- like structures ^, ^, are 
still displayed by their random counterpart, which 
are therefore good candidates to model quasicrys- 
talline materials 0. The best description for real 
quasicrystals remains an open question. 

Two years ago, we presented initial results on two- 
dimensional random tilings of large codimension m 
in which, as had been anticipated by Henley 
an analytic approach is more likely to be developed 
than in the finite-dimensional case. The present pa- 
per contains preliminary results in the understanding 
of these systems. We develop a mean-field theory, 
based upon an iterative walk-on-tiling construction 
of large codimension tilings, providing a new insight 
into this infinite codimension limit. Except for the 
diagonalization of an operator, all this work relies on 
exact calculations. The results so-obtained are val- 
idated by Monte Carlo simulations, based upon the 
same iterative process. A more complete and explicit 
presentation will be published elsewhere ^ . 

In the two following sections, we briefiy define 
what infinite codimension tilings are and we recall 
results from reference [Q. In section 4, we explain 
the mean-field theory and we compare it to results 
of Monte-Carlo simulations. 
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2. Model 

We define two-dimensional random tilings of rhombi, 
denoted as L* ^ 2 tilings, as projections from a 
£)-dimensional cubic lattice into a two-dimensional 
plane. Facets of the cubic lattice project into 
rhombi [|[ |l^, |ll|. The projection is chosen so that 
the rhombi have unit edge length and their angles 
are multiple of -k/D. As a consequence, the tiles 
are arranged in tilings with statistical 2i?-fold rota- 
tional symmetry. Figure ^ displays such tilings. In 
this case, they fill regular 2£'-gons. Such polygonal 
boundary conditions will be called fixed in the fol- 
lowing. We will also consider free boundary tilings 
in this paper. In the de Bruijn grid representa- 
tion H 
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|14| ], any two lines of different orienta- 
tion cross in the fixed boundary case, whereas they 
do not necessarily cross in the free boundary case. 




Fig. 1. Two examples of fixed boundary tilings for D — 5 
and 7. In this case, the side lengths of the 2_D-gons are 
all equal to p = 8. 

Before going on, let us mention that the effect of 
boundaries was recently studied and turned out 
to be highly non trivial: fixed boundaries have mi- 
croscopic effects on tilings at the infinite size limit, 
which are characterized by a gradient of entropy be- 
tween the boundary and the center of the tilings. As 
a consequence, fixed boundary tilings have a lower 
entropy than free boundary ones. Note that periodic 
boundary conditions give the same entropy as free 
ones at the infinite size limit fl^ , and that this en- 
tropy only depends on the relative fractions of tiles 
and not on the sample shape. 

3. Thermodynamic limit 

We demonstrated previously 0] that when the codi- 
mension|^ of the plane tilings becomes large, the en- 
tropy difference between fixed and free boundary 
conditions vanishes. More precisely, the structural 
inhomogeneity between regions close to the bound- 
ary and the bulk observed in finite codimcnsion 

^That is the difference D — 2 between the cubic lattice 
dimension D and the tiling dimension 2. 



tilings tends to disappear. Moreover, it was argued 
that this limiting entropy is still valid when the side 
length p of the polygonal boundary remains finite as 
the codimension goes to infinity. In particular, one 
can choose p = I m order to perform more simple 
calculations [0. In this reference, an upper bound of 
the limit entropy, a < log 2, was also given. 

4. Mean-field theory 

Both the mean-field theory and the numerical calcu- 
lations presented below rely on an iterative construc- 
tion of large codimension tilings, based upon the de 
Bruijn representation of tilings [^2j. This point of 
view can be used either in the fixed boundary con- 
text j^], or in the free boundary one. 

Walk-on-tiling algorithm 

As illustrated in figure ||, a 13 — > 2 tiling can be 
seen as a collection of non-intersecting paths chosen 
on a 13 — 1 2 random tiling. These paths are 
directed from bottom to top. They follow the edges 
of the D — underlying tiling. Then these paths 

are "opened" along a new edge direction in order to 
form the D-th de Bruijn family and then a Z3 ^ 2 
tiling. Conversely, de Bruijn lines of family D in 
a, D 2 tiling can be collapsed into paths on a 
D — 1 2 tiling. As a consequence, there is a one- 
to-one correspondence between 13 — > 2 tilings with 
p de Bruijn lines in family D and collections of p 
non-intersecting directed paths on _D — 1 — ^ 2 tilings. 





Fig. 2. Iterative process for the construction of D — + 2 
tilings: paths are chosen on a D — 1 ^ 2 tiling (left), and 
are "opened" in a new direction to form de Bruijn lines 
of the D-th family (right). 



Mean-field theory and Monte Carlo simula- 
tions 

As we explained above, as D becomes large, the 
enumeration problem can be simplified into a prob- 
lem with a single line in each de Bruijn family. In 
other words, de Bruijn lines of the same family do not 
interact, and the number of p-line configurations can 
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be factorized into a product of independent single- 
line configurations. As a consequence, the entropy is 
equal to 

. . lim lim ^HiiM^, (1) 

where iV£)(fc) is the number of /c-steps paths on a 
Z3-dimensional random tiling || . 

Now, the mean-field approximation assumes that 
the steps of such paths are uncorrelated. Then the 
number of paths reads: 

fc 

Nnik)^l[c„ (2) 

where the Cj are the numbers of choices to be made 
at each vertex j to add a new step to the path (see 
figure H) . 




Fig. 3. A random path on a tiling as a succession of 
uncorrelated vertices (the bottom-to-top path has been 
rotated into a left-to-right path for sake of convenience). 

Keep in mind that, even if they are uncorrelated, 
the vertices j belong to a random £) — > 2 tiling. 
Therefore the numbers of choices Cj are distributed 
according to a probability distribution TToic). When 
D tends to infinity, this distribution tends toward a 
limiting distribution, denoted by 7r(c). Thus, after a 
short calculation, equations (0) and (§) become: 

oo 

= ^'^ic)\ogc. (3) 

c=l 




Fig. 4. A path going through a vertex s of a D — > 2 tiling 
gives birth to two vertices si and S2 of a D + 1 ^ 2 tiling. 

To get the mean-field distribution 7r(c), let us first 
denote by tt^i {q, p) the fraction of vertices on a _D ^ 2 
random tiling with q in-coming edges ( "legs" ) and p 



out-going ones ("arms") (see figure ||). The distri- 
butions 7r£)(c) and 7r(c) will easily be derived from 
■K£){q^p). If no correlations are taken into account, 
a (q,p)-vertex s of a D — > 2 tiling will be visited 
by the paths with a probability 7Td{QtP) and will 
then give birth to two vertices Si and S2, as illus- 
trated in figure IJ. These latter vertices belong to a 
_D 4- 1 2 tiling. Each leg and each arm of s will be 
chosen with probabilities 1/q and 1/p, respectively. 
As a consequence, in this mean-field approximation, 
the probabilities ttd+i{<1tP) can be written as linear 
combinations of the probabilities T:D{q,p)- The cor- 
responding linear operator will be denoted by A. It is 
infinite-dimensional since q and p can be arbitrarily 
large when D goes to infinity. 

The limiting distribution T:{q^p) is the fixed point 
of A, that is the eigenvector associated with the 
eigenvalue 1. We have calculated this fixed point 
numerically, yielding the values of 7r(c) listed in ta- 
ble |l|. The mean-field entropy (^) is then approxi- 



c 


1 


2 


3 


4 


Mean-field 


0.30 


0.45 


0.19 


0.04 


Numerical 


0.26 


0.51 


0.20 


0.03 



Table 1. The first values of the limit distribution n{c), 
obtained both in the mean-field approximation and nu- 
merically, by Monte Carlo simulations. 

matcly equal to cr = 0.60. The distribution ■n{q,p) 
also provides the mean-field correlations between p's 
and g's: their covariance is cov{q,p) ~ — 0.44(^. 

Let us mention that we have not taken into ac- 
count the fact that at the step D of the iterative 
process, only a fraction of order 1/D of the vertices 
are visited by paths This modification does not 
alter the fixed point of the process but only the 
finite D corrections: they have a power-law behavior 
instead of an exponential one. 

To close this section, let us emphasize that these 
mean-field values are satisfyingly close to numerical 
values obtained by Monte Carlo simulations as 
displayed in table |l|. The corresponding numerical 
entropy and correlations are approximately equal to 
0.57 and -0.36, respectively. 

5. Conclusion 

The above mean-field theory provides quite satisfy- 
ing numerical results concerning the local structure 
of large codimension random tilings, in terms of ver- 
tex statistics distributions. As a consequence, it pro- 
vides a good approximate value of the limit entropy 
of such systems. 

^As a matter of fact, this value seems to be equal to -4/9. 
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To go further and to get better approximations, it 
will be necessary to include the role of spatial cor- 
relations in our calculations, by taking into account 
the distribution of tiling patches bigger than single 
vertices. 

Moreover, it would be useful to get informations 
about finite D corrections to the limit entropy cr, 
since finite codimensional systems are related to real 
quasicrystals. Two ingredients will be taken into ac- 
count: the power-law corrections of the mean- field 
theory as discussed above, and the effects of contacts 
between dc Bruijn lines, which tend to decrease the 
entropy 0. 
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